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C3 ' Abstract 

^3 : 

1-^ , We propose a simple extension of the minimal supersymmetric standard model 

Q^ ' by introducing a gauge singlet in addition to right-handed neutrinos. The model 

■ resolves the strong CP problem by Pecci-Quinn symmetry, explains the origin of 
,— i! left-handed neutrino masses as well as MSSM /^-parameter. It also gives rise to 

I thermal inflation, baryogenesis and dark matter in a remarkably consistent way. 

I ■ Interestingly, resolution of moduli problem by thermal inflation constrains tightly 

axion coupling constant and flaton decay temperature to be fa ~ 10^^ GeV and 
r~| \ T(j ~ 100 MeV, respectively. Model parameters in this case are likely to give right 

amount of baryon asymmetry and dark matter at present. The main component 
^SJ . of dark matter is expected to be the axino whose mass is nearly fixed to be about 

> : 1 GeV. 

^ : 

(N . 

(N ; 1 Introduction 

^' ■ 

' Although its great success in describing particle physics, the standard model (SM) 

■ of particle physics has been faced on some big questions: hierarchy problem associ- 
• • , ated with the unnatural stability of weak scale Higgs mass against large quantum 

. ^ \ corrections induced by the large hierarchy between electroweak (mew ~ 10^ GeV) 

^ ■ and Planck (Mpi ~ 10^^ GeV) scales, strong CP problem [ ] and the origin of 

left-handed neutrino masses which is implied by the neutrino oscillation [-:]. Any 
plausible theory for particle physics, which is based on SM, should address all 
these questions. 

The minimal supersymmetric standard model (MSSM) is a very natural and 
simple supersymmetric extension of standard model, in which large hierarchy prob- 
lem is absent by the help of supersymmetry (SUSY) [ ], hence it is quite attractive. 
MSSM may be the right direction to a theory for low energy physics, but it has a 
mysterious parameter /i, the dimensionful coefficient of Higgs bilinear superpoten- 
tial term, whose origin need to be explained [ 1]. In addition, MSSM still should 
be extended to resolve strong CP problem and explain the origin of left-handed 
neutrino masses. 

Meanwhile, string/M theory, which may be the fundamental theory at high 
energy, predicts the existence of light particles with gravitationally suppressed 
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interactions: gravitinos and moduli/modulinos. Typically, those particles are ex- 
pected to be produced too much due to large reheating temperature or coherent 
oscillation after primordial inflation, disturbing the successful Big-Bang nucle- 
osynthesis (BBN) or over-closing the present universe ["), (>]. The most compelling 
solution to this gravitino/moduli problem is thermal inflation [ ] which is very 
likely to occur in the framework of SUSY. But thermal inflation invalidates most 
of known baryogenesis mechanisms, which work typically above or at the elec- 
troweak scale, by diluting out pre-existing baryon/lepton asymmetry ^. Since 
the reheating temperature after thermal inflation is typically well below the elec- 
troweak scale, baryogenesis mechanism is difficult to work in general. 

In this paper, we propose a simple extension of MSSM by introducing just a 
gauge singlet and right-handed neutrino chiral superfields, but without any ad- hoc 
mass parameter except soft SUSY-breaking ones. In the framework of gravity- 
mediated SUSY-breaking, the model realizes Peccei-Quinn symmetry [10, f 1, 12] 
to resolve strong CP problem and explains the origin of left-handed neutrino 
masses through see-saw mechanism [13] as well as that of MSSM /z-parameter in 
a very efficient way. The model also provides remarkably consistent cosmology 
including thermal inflation, baryogenesis and dark matter. 

This paper is organized as follows. In Section 2, we propose a model and de- 
scribe how it provides proper low energy physics. In Section 3, we briefly describe 
cosmology which the model can realize. In Section 4, we conclude. 

2 The model 

Motivated by the drawbacks of MSSM, strong CP problem, the origin of left- 
handed neutrino masses and fj, parameter, we extend the MSSM to the following, 

assuming SUSY-breaking is mediated by gravitationally suppressed interaction ^ 

3 . 

W = XuQHuU + XdQHdd + XeLHae + X.LHuU + A^^^^^ + ]-X^^u^ (1) 



^One may think that AfHeck-Dine baryogenesis in a scenario of gauge- mediated SUSY-breaking can 
work well even in the presence of thermal inflation [S]. However, the formation of Q-balls makes it 
difficult to work [ ]. 

^In the scenario of gravity-mediated SUSY-breaking, the symmetry breaking scale of thermal in- 
flation required to resolve moduli problem coincides accidentally with that of Peccei-Quinn symmetry. 
Hence, in the scope of this paper, it is natural to consider such a mediation scenario as the framework 
of our argument. We refer the reader to ref. [ '] for the case of mirage- mediation of SUSY-breaking 
[15], and will discuss the cases of other mediation scenarios, for example gauge-mediation [IG, 17], in 
other place. 

^For the MSSM /i-term, we can use a renormalizable interaction A^$i?„i?d instead of the non- 
renormalizable interaction Xfi^^HuHd/Mp\ which we are using in this paper. Our model then involves 
only renormalizable interactions though ^ 10^^ is hierarchically small. It may be still plausible to 
consider such a hierarchically small Yukawa coupling, since Yukawa couplings of MSSM already have 
a hierarchy. Compared to the case of non-renormalizable interaction for /i-term, the physics is barely 
affected and can be applied to most of known mediation mechanisms of SUSY-breaking except gauge 
mediation. 
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with an assumption 



where and are respectively the soft mass-squared parameters of L and 
around electroweak scale. In Eq. (1), indices for gauge group and family structure 
has been omitted for simplicity, Mpi = 2.4 x 10^^ GeV is the reduced Planck mass, 
<I> = + \f2Qa + • • • is a gauge singlet chiral superfield and v is the right-handed 
neutrino superfield. Various new Yukawa couplings are constrained by low-energy 
physics which will be addressed shortly. Note that the model has an accidental 
{7(1) symmetry. Normalized such that ^ has charge one, the charges assigned to 
fields are as follows. 

C/(l) {g, L, Frf, n, d, e, v] = {\, 3/2, -1, -1,0, 0, -1/2, 1, -1/2} (3) 

The key feature of our model is that if A$ ~ 0(1) the associated Yukawa inter- 
action can drive the soft mass-squared parameters of </» and right-handed sneutrino 
(denoted as v except cases of confusion) to be negative through renormalization 
group running ^. If v develops vacuum expectation value (vev) before ^ does, our 
model becomes inconsistent with low energy phenomenology since MSSM ^-term 
is not reproduced. In order to avoid this disaster, we assume at Planck scale 



m 



(4) 



where and are respectively the soft mass-squared parameters of (/> and 
and is the ^-parameter of the trilinear coupling associated with A$. In this 
case, only ^ can develop non-zero vev and the spontaneously broken accidental 
{7(1) symmetry can be identified as the Peccei-Quinn (PQ) symmetry for DFSZ 
axion [12]. 

The potential along is of the form 



F(|</.|):^yo + m2(Q=|A$0|)|0p 



(5) 



where m|((5) is the running soft mass-squared of </> evaluated at a renormalization 
scale Q. Now 



dV 



2 + 
1 + 



d 



dln[ 
d 

dlnl. 



^(101) 
2 + 



din 



(6) 
(7) 



and the renormalization group equation of is 
dmi I 



dlnQ 87r2 
hence at vacuum where dV/d\(j)\ = we find 



-^\\^mml + ml + \A^\^) 



(8) 



m,^(<^o 



167r2 



|A$p {ml + |^$p) 



(9) 



*See Appendix for renormalization group equations of relevant parameters. 
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where the parameters at right-hand side are evaluated at Q = \X^(j)o\ and cpo is the 
vacuum position of (p. Since the renormalization group runnings of X^, mf, and 
are rather slow as shown in Appendix, from Eq. (8) we may approximate 
m3,(0o) crudely as 



mliM ~ m^(Mpi) + ^|A^p {ml + \A^\') In 



(10) 



Then the vacuum value (j)o, which is assumed to be real without loss of generality, 
is 

/ifo, r 1 11 

(11) 



■ exp 



1 1 

2~ a 



where 



A<j) 

= 87r2 m2(Mpi) 
and requiring zero cosmological constant, one finds 

Vo = ^am2(Mpi)<^2 



a 



(12) 



(13) 



Decomposing (j) as (p = [(po + s / \/2^ exp [ia/ (\/2(/>o)], the physical mass of 
saxion (s), denoted as mpg, is 



1 fv_ 

2 



dm\ 



00 



din I 



00 



1 

8^ 



\X^\^ml + \A^\^) 



(14) 



The masse of axon (a) at zero temperature is [18] 

, , /l0i2GeV\ 
ma ~ 6 X 10~^ eV \ j J 

where the axion coupling constant fa is defined as 

\/2(/)o 



fa 



N 



(15) 



(16) 



with = 6 the coefficient of U (I)pq-QCD anomaly in our model in normalization 
of unit PQ charge for PQ field (0). fa is lower-bounded by the cooling rate 
of SN 1987A [ ] and upper-bounded by over-closure limit of cold axions from 
misalignment and strings [18]. If there is no dilution after axion condensation 
from misalignment is formed, currently allowed window of fa is given by 

10^ GeV < < 10^2 GeV (17) 

The masse of axino (a) is generated at 1-loop in our model. It is given by [20] 

1 



rrin 



167r' 



■A|A$ ~ 0.6 GeV A| 



100 GeV / 



(18) 



hence axino is the lightest super symmetric particle(LSP) in our model. 
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The large vev of (f), implied in Eq. (17), also makes v become very heavy due 
to the Yukawa coupling associated with A$ in Eq. (1), and integrating out v in 
Eq. (1) gives an effective low-energy superpotential which contains a left-handed 
neutrino mass term: 

T^eff = KQHuU + XdQH,d+ X,LH,e + X^^^^ - J^^r^^ (19) 

Mpi 2 X^cpo 

In order to provide right size of masses to left-handed neutrinos, we need 

v2„,2,:„2fl / \ \2/in-2„Ar\ / 1 nl2 



A, = . 0.33 (^] ('^) C^^) sin'p 



(20) 



where vfiy^ is the left-handed neutrino mass, v = 174 GeV is the vev of Higgs 
field and sin/3 = Vu/v with Vu the vev of up-type Higgs field (Hu)- Note that in 
Eq. (19) MSSM n parameter is given by 

Note also that using Eq. (20), we can re-express Eq. (11) as 

"^^(^P') ^ /_l_l^Mo (22) 

+ |^$|2 2V2^V 2 Mpi 

/ \ VlO-2eV\ /10i2GeV\ , , 

~ 0.14 ■ sin^ /3 (23 

which provides a constraint on mass parameters. 

As described before this, introduction of a gauge singlet chiral superfield allows 
natural realization of the PQ-symmetry to resolve the strong CP problem, and 
explains the origin of small mass of left-handed neutrino and MSSM //-parameter. 

Our model is similar to the models of refs. [ ] studied in refs. [2l ] and 
those of refs. [23, 24], but it differs from those models mainly by how symmetry 
breaking field is stabilized and/or how a right-handed neutrino obtains large mass 
for a see-saw mechanism. In our model, instead of non-renormalizable higher 
order operator, the effect of field-dependent running of a mass-squared parameter 
is used to stabilize the symmetry breaking field whose vacuum expectation value 
provides a large mass to right-handed neutrinos. Although it may be thought of 
a rather simple variation of the models, the difference makes our model simpler 
and maybe more natural with rather different cosmology which will be described 
in the next section. 



3 Cosmology 

As our ansatz, we assume only HuHd fiat directions [_ >] and (p develop 

non-zero field values temporarily or eventually ^. This will be justified shortly in 
subsequent arguments. Parametrizing LHu and HuHd fiat directions as 

L={Q,lf , Hu = {hu,^f , Hd = {Q,hdf (24) 

^Although HuHd may be stable near the origin as we are assuming here, dynamics of fields can lead 
non-zero field value of the flat direction. 
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with the remaining I?-term constraint 

D = \K\^ -\hd\^ 
we find a potential from Eqs. (1) and (5) 



V = ml\lf 



huhd + c.c. 



A, 



K 



1 



+ 



\hr. 



^Mpi 



+ 



2A„ 



■Mpi 



(25) 

(26) 
(27) 

(28) 
(29) 



where terms in the first to forth lines are soft SUSY-breaking mass terms, A-terms, 
F-terms and D-terms, respectively, and = {gl +g2)/^- Based on this potential, 
in this section we will describe cosmology of our model including thermal inflation 
which may be indispensable to resolve moduli problem, baryogenesis and dark 
matter. 



3.1 Moduli problem and thermal inflation 

In string/M theories, that might be candidates of the fundamental theory, there are 
flat directions called moduli. A modulus {(p) has interactions suppressed by Planck 
scale, and starts to oscillate coherently when expansion rate becomes comparable 
to the mass of the modulus, rriip. We assume that radiation is dominant at this 
time, then the abundance of the moduli (called Big-Bang moduli) is 

!V) .(^^n)"^45xlO'fiISX)"= (30) 

s Jbb \mp J \ J 

where is the number density of moduli and s is the entropy density. The decay 
rate of moduli is given by 

where 7^ ~ C'(l) is a numerical coefficient, hence moduli decay after BBN. In order 
for successful BBN, the abundance of moduli when they decay is constrained to 
be [26] 

^ < l°-"G<^V , 10-" f ilSX) (32) 

s nip \ "^ip J 

Therefore, we need a dilution more than 0(10^*^). 

Thermal inflation is the most compelling solution to this moduli problem. 
Moduli may start to oscillate coherently while (p is held around the origin due to 
the interaction with v which is in very hot thermal bath. Thermal inflation begins 
when the potential energy density at the origin becomes dominant over the energy 
density of moduli at a temperature given by 



6 



where we have used ~ m^Mpj as the initial energy density of moduh coherent 
oscihation. Thermal inflation ends as (p starts to roll out when temperature drops 
to the critical temperature 



In 



"T-soft 



1/2 



(34) 



where mgoft ~ 1 TeV is the scale of soft SUSY-breaking. 

The coherent oscillation of (p after thermal inflation eventually decays with 
rate 

r,/, ~ T^^aa + r<^^SM (35) 

where 



1 m 



PQ 



647r 



(36) 



is the partial decay rates of cj) to axions. r^_i.gM is the partial decay rates of (p to 
SM particles. It is mainly due to the mixing between cp and Higgs fields, induced 
by the term A^^^^^ifrf/Mpi. Since mpq ~ 0(10 - 100) GeV for A$ ~ 0(1) from 
Eq. (14), r0_^sM is expected to be dominated by the decay to bottom quarks and 
given by [24] 



i-SM 



1 

in 



\B\ 



m , 



m 



f\ 



mt 



m 



(37) 



where tjia is the mass of CP-odd neutral Higgs particle, i? = ^„ in our model. 



fix) 



ex 



(1 



y) 



(38) 



and 



12m? 
e ^ ~ 0.02 

mf 



(39) 



with m-fe and m/j the masses of bottom quark and light neutral Higgs particle, 
respectively. Successful BBN limits energy contributions of relativistic non-SM 
particles, so it requires T^^aa/^(j)^SM ^ 0.3 [24] which is easily satisfied for mpq ~ 
0(10 - 100) GeV. From Eqs. (36) and (37), the decay temperature of flaton {(p) 
is given by 



g*{Td 



26 GeV 



-1/4 



?vr4g,(rd 



(r</>-i>sMr, 
1 - —k- 



,1/4^1/2 



PI 



\B\ 



m , 



10^2 GeV 



mpQ0o 



rat 



m 



I V ITeV 



PQ 

2 



30 GeV Y 

"^PQ / 



ml 



m 



PQ, 



(40) 
(41) 



;42) 



where we have used g*{Td) 
of mpq/m/i. 



100 in the last line. Fig. 1 shows function 
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10' 
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10" 




mpQ / rrih 



Figure 1: Flaton decay temperature versus flaton mass: Td versus mpg/m/i for 0o 
10^2 GeV, = 10^ GeV, ruh = 125 GeV and uia = 2\B\. 



The e- foldings of thermal inflation is, from Eqs. (33) and (34), 
= In^ 

c 

4 



~ 7.6 + ^ In 

-im 

2 



30 GeV y / lTeV \ 
1012 GeV y V mpQ ) \ ) 



ITeV/ 



_10i2 GeV, 

and the dilution factor due to the entropy release in the decay of (p is 
Vo 



A. 



r3Td 

2 X 10^^ 



1012 GeV 



\V30GeV\ /10-iGeV\ , / 



(43) 



' lO'^msoft \ 



(44) 



(45) 

-3/2 

(46) 



|A<i.(/>o| / 

This is little bit small to dilute out moduli produced before thermal inflation to a 
safe level. 

Fortunately, there is additional source of dilution. We assume 



Tc < Tlh,. 



(47) 



with 



Then, as temperature drops below TiHu-> which is expected to be held near 
the origin at high temperature rolls out from the origin before flaton is destabilized. 
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It is stabilized by the F-term potential, {X^lhul"^ in Eq. (28) at 

^ \hu*\ ~ mLH„/\K\ (49) 

The initial energy density of the coherent oscillation of LHu is 

Vlh^ ~ mi^J/,|2 ~ mi^JlKl^ (50) 

If this energy can be quickly dumped into radiation we would get an additional 
dilution to Eq. (45). 

The condensation of LH^ decays initially through prheating [_ ; ] and eventu- 
ally through perturbative decays. For X,^ > 10^^, the perturbative decay of the 
condensation is dominated by the coupling Aj/ with a partial decay rate given by 

6 

^LH^ ~ ■^\Xu\'^rnLH^ (51) 



Using Eqs. (13), (12) and (20), one finds 
Tlh^ 280 



^Tti sin2/3 VlO-2eV7 V 

where ffxi is the Hubble expansion rate during thermal inflation, hence the co- 
herent oscillation of LH^ decays within a Hubble time. 

It is expected that preheating is at best dominant for modes k ~ miHu and 
order one fractional energy density of LHu is still in the form of condensation, 
hence the Universe is slightly reheated to a temperature 

T*-'Vli,l^\X,\-'/^mLH^ (53) 

before flaton is destabilized. Thermal inflation is then extended by e-foldings given 
by 

Nlh^= In ^ ^ln|A,| (54) 

and, using Eq. (53), the dilution factor due to the entropy release in the decay of 
LHu is 

y3/4 

Alh^ \X.\-^/^ (55) 

Therefore, from Eqs. (30), (45), (55) and (20) the abundance of Big-Bang moduli 
when they decay is 

1 3/2 



10^2 GeV^ ^ 



) V30GeV/ VlO"^GeVy ^\\X^(t)Q\) 



3/2 

(57) 



^If a coupling of LH^ to MSSM particles is larger than Aj/, the particles becomes heavier than LHu, 
hence the decay of LH^ is kinematically forbidden. 
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which can manage to satisfy the bound from observation in Eq. (32). Thus, Big- 
Bang moduh can be diluted out to a safe level. 

Moduli are also produced after thermal inflation due to the finite energy density 
of thermal inflation. The abundance when they are produced is 



T/2 



s Jti T^m^M^^ 
Using Eqs. (45) and (13), the abundance when they decay is 
(n^/s)Tj VoTd 



(58) 



(59) 



Compared to Big-Bang moduli of Eq. (57), it is negligible, hence it is harmless as 
long as Big-Bang moduli is diluted to a safe level. 

Apart from the dilution of unwanted relics, thermal inflation has a very in- 
teresting aspect. It wipes out pre-existing gravitational wave backgrounds and 
generates its own one which may be detected at BBO or DECIGO type experi- 
ment [2S]. Since the peak frequency and amplitude depend on the details of model 
parameters and physics of phase transition, careful investigation may be necessary 
to see if the wave in our model is really detectable. 



3.2 Baryogenesis 

As (/) rolls out below Tc, thermal inflation ends and v becomes very heavy. By 
integrating out v in Eq. (1), one finds an effective superpotential 



W = XuQHuU + XdQHdd + XeLHde + A^— — - 

Mpi 2 A^(p 

which gives a potential 



V = ml\l\^ + mljK\^ + mjjjhdl^ + ml 



b^Khd 1 , XI {IKf 



Mpi 2 
4>^hd XlPhu 



^ Ml 



PI 



+ 



A<j>i;^ 
(l?hu 



+ c.c. 



PI 



+ 



^A„, — — h 



^ Mpi 2 A$(/)2 



2\2 



(61) 

(62) 
(63) 

(64) 
(65) 



for (p > msoft- 

As (p becomes large, LH^ is shifted to a larger value 



\l\^^\hu\ 



\X^^\ i^^l2mlH^ + \A^\^ + \A^ 



(66) 
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For msoft <C (/> <C (^i'Oi the phase of LH^ is determined by the Ai^-term in Eq. (63). 
As (p reaches its vev, MSSM ^-parameter is reproduced, hence LHu is Hfted up 
and brought back into the origin ^. Simultaneously, HuHd becomes large ^ due 
to the cross term of the third term in Eq. (64), and the cross term of the first 
term in the same equation provides angular kick to LHu- This is nothing but an 
Affleck-Dine (AD) mechanism [29] to generate charge asymmetry. 

The generated lepton asymmetry is expected to be conserved due to the rapid 
preheating of the AD fields [ ]. The oscillating AD fields decay earlier than 
flaton, reheat the Universe partially and activate sphaleron process [ ] in which 
lepton asymmetry is converted to baryon asymmetry. Eventually, flaton decays 
and reheats the Universe in order for successful BBN. There is huge amount of 
entropy release in the flaton decay. The baryon asymmetry at present is then 
given by 



/I I \ 2 

UB riB Td riL riAD riL rriLHu ( 



mpQ HAD uipq had rnpQ \(po / mpq 



(67) 



where n^, ni and had are number densities of (p, lepton asymmetry and AD fleld 
respectively, and Iq is the value of / when (p reaches (po from the origin. From the 
experience of a similar model of Ref. [24], we expect 



"AD 

and from Eqs. (66) and (20) 



hence we expect 

riB _ -^Q-io f riL/n AD 



s K 10-3 

X (H^Y /SOGeVy / 10-2eV \ /lO^^GeVV / \ 
VlTeV; \ mpQ J V 7 V '^0 / VlOOMeV/ ' 

Note that the model parameters to resolve moduli problem give right order of the 
baryon asymmetry which can match the present observation. 



3.3 Dark matter 

In our model, axinos and axions are expected to be the main components of dark 
matter at present. These dark matter components can be cold, warm or even hot, 
depending on their masses and how they are produced. In this subsection, we will 
derive cosmological constraints from the dark matter. 

^In order to avoid possible local minimum along LHu, we need mj^+mj^ + ImP > l^i/ otherwise 
LHu may be trapped there that is disastrous. 

^ Large HuHd holds dangerous quark and lepton flat directions near the origin, hence considering 
only LHu, HuHd and (p is justified. 
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Axinos produced in the flaton decay As the LSP, axinos are produced 
mainly in the decay of flaton with a rate given by [24] 



327r(/>2 



(71) 



where Ua = {dlnrria) / (dln|0|). In our model, from Eqs. (18), (104) and (105) 

5 



On 



87r2 



The axino number density from the flaton decay is 



2T, 



Ha 



a{t'fp^dt' 



(72) 



(73) 



(74) 



where a is the scale factor and p^f, is the energy density of (j). Using Eqs. (40) and 
(81), one find [24] 

l^d 2.2T(jr0^aa 

s m^TsM 
From Eqs. (40) and (71), the current axino abundance is 

,8 f \ n-d 



^d 



5.6 X lO'' 



GeV/ 



(75) 



.1/2 



~ 0.36- 



10 



.1/2 I ^1/2.^ s 
SM \y* V d/ 



an 



10- 



100MeV \ / lO^^ GeV \ 
IGeV; V^^J ) V ^0 ) 



2 

76) 



Therefore < f^CDM — 0.25 requires 



m-a < 0.9GeV 



rl/2 I 



10 



100 MeV 



10" 



1\ 2 



an 



1012 GeV 



(77) 



which can be easily satisfied for A$ < 1 from Eq. (77). 



Axinos in the decay of thermally generated NLSP If the next-to- 
lightest supersymmetric particle (NLSP) denoted as x is of neutralino type, the 
coupling X^cfP' HuHd/ M-pi in Eq. (1) leads to the decay of x to axino via various 
channels with rate [31] 



167r ' 



(78) 



where Ca ~ 1 may contain a factor of / and we have neglected the masses 
of decay products. 



The thermal bath generates x with the number density 

= ^ I dk — 



1 /"^ „ k'^ 



-7r2 , 
^ ./O 



exp v/ + 1 



(79) 
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and they subsequently decay into axinos. The late time axino abundance is thus 
estimated as 



= ^ / a{t'fn^{t')dt' (80) 



whose numerical solution results in [24] 

1 /2 

— = 9* {Td)g* iTx)[— 1 Pa{x) (81) 



where Fa{x) is approximated as 



with 



5.3x'' for a; ^ 1 , , 

5.4 for X > 1 ^ ' 



3 \a*{Tx)j 

and Ty. ~ 2m^/21 being the temperature at which the axino production rate is 
maximized. From (78) and (81), the current abundance of axino dark matter is 
obtained 

2 7C-r 1 i 1 AO^^GeV y 

y^fT.lo^/^my VlO^GevJllGevJ V 'Ao ) ^ 

(85) 



where we have used ~ F(^^sM- Therefore, for x <C 1, one finds 
which does not exceed ilcDM if the flaton decay temperature satisfies 



< 7.5 GeV 



100 GeV/ V9*(7d 

1/4 .^N 5/4^ .-.„2 



^-i/90ZiiM_^ /10^\ / lGeV \ / cPo Y 

V 103 H / V y vio^^Gevy 



1/7 

(87) 



Note that in order to resolve moduli problem (See Eq. (57)) and explain baryon 
asymmetry of the Universe (See Eq. (67)) we need ~ 0(100) MeV, and in that 
case Eq. (87) is satisfied automatically. Since is expected to be much less than 
the freeze-out temperature of which is about m^/20 [ ], in order to avoid 
direct production of x from the flaton decay, we may require < 2my, which 
seems to be typical in our model. 

Axinos will also be produced by the decay of x after they freeze out. However, 
the standard Big-Bang neutralino freeze-out abundance is good match to the dark 
matter abundance, our freeze-out abundance of x will typically be less than the 
standard abundance, and nia ^ m^, therefore the axino abundance generated 
after the freeze-out should be safe. 
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Axions from misalignment and strings The current abundance of axions 
produced in these ways is given by [ ] 



" Aa Vl0i2GeV 

where is the dilution factor of axion misahgnment due to the late time entropy 
release of thermal inflation. For <^ 1.3 GeV [24], 

1.3 GeV^ '-^^ 



A, ~ ( ) (89) 



Td<2.0GeV(^^^^) (90) 



hence for fl^^'^ < J^cDM — 0.25 we need 

0^2 GeV ^ 

Since we expect < O(100)MeV, these axions are likely to be a sub-dominant 
contribution to cold dark matter at present. 

Axions in the flaton decay Axion can be warm or hot when they are pro- 
duced in the decay of flaton. The constraint on hot dark matter comes from 
CMBR and structure formation [32, 33]. Currently allowed hot dark matter frac- 
tional contribution to the present critical density is ^hdm ^ 10"^ [18]. The 
constraint may be more stringent as suggested from the analysis of the early re- 
ionization of the Universe at high redshift [ ]. Taking into account the recent 
analysis of WMAP 5-year data [•> <], the allowed warm/hot dark matter fractional 
contribution to the present critical density is likely to be 

^^WHDM < 10-3. (91) 

We will take this as the upper bound on the fractional energy density of our 
warm/hot dark matter. 

Axions produced by the flaton decay have a current momentum 

= (!.2) 

ao 2 

where a is the scale factor at the time they were created and oq is the scale factor 
now. Whereas, the current momentum of an axion produced at td = is 

Pd = ^ (93) 

ao 2 
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where and 5f are respectively the total entropy at decay time td and present, 
so it would be highly relativistic now. The current number density spectrum is 
given by 



dpa \aoJ H ml Hp\ 



which may provide an observational test of our model in the future. The energy 
density of the axions is 

Therefore, assuming that the hot axions are still relativistic now, their current 
energy density is estimated as^ 

^ 4.3 X 10-5 1^ (98) 

\g*\-l-A)J l(/)-i^SM 



4 Conclusion 

In this paper, we have proposed maybe the simplest extension of MSSM by intro- 
ducing just one additional gauge singlet chiral superfield and right-handed neu- 
trino superfield {v) to MSSM without any ad-hoc mass parameter except soft 
SUSY-breaking ones: 

W = XuQHuU + XdQHdd + XeLHde + KLHuV + A^^^^^^ + ^A^^i/^ (gg) 

ikZp] z 

with an assumption for the soft mas-squared parameters of L and Hu 

ml + m\^<^ (100) 

and surely ml + m?^^ + > around electroweak scale. This model can realizes 
Peccei-Quinn symmetry to resolve strong CP problem and explains the origin 
of left-handed neutrino mass and /^-parameter of MSSM. In addition, the model 
realizes thermal inflation followed by Affleck-Dine type leptogenesis and provides 
dark matter which matches well to observation in a remarkably consistent way. 

At low energy, the scalar component of the gauge singlet field, (j) has negative 
mass-squared around the origin due to strong Yukawa coupling to hence devel- 
ops non-zero vacuum expectation value. In the absence of self-interactions, the 
field can be stabilized radiatively around intermediate scale, so it can be identified 
as the Pecci-Quinn field which breaks the Pecci-Quinn symmetry spontaneously. 
The vacuum expectation value of (j) make v heavy, hence integrating out v at low 
energy generates the mass term of left-handed neutrino through seesaw mecha- 
nism. The vev of (f) also generates MSSM /i-parameter. 

^The energy density of thermally produced axions is i^a ^ rria/lSleV, hence it will be subdominant 
[30]. 
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Cosmologically, this model gives rise to thermal inflation while (j) is held near 
the origin due to large thermal effect at high temperature. Notably, the non-trivial 
assumption of m\ + w-^^ < allows LHu flat direction to be destabilized before 
the end of thermal inflation. This causes an extension of thermal inflation by few 
e-foldings. Although it is small, the additional e-foldings are crucial to dilute out 
moduli to a safe level. Thermal inflation ends as is destabilized at a critical 
temperature. As </> reaches its vacuum expectation value, MSSM //-parameter 
is reproduced, leading Affleck-Dine leptogenesis involving LHu and HuHd flat 
directions. The oscillation of Affleck-Dine fields is expected to be quickly damped 
due to a rapid preheating induced by oscillation and preheating of flaton field (j), 
so the generated lepton asymmetry can be conserved. 

The eventual decay of flaton reheats the Universe with a temperature ~ 
O(IOO) MeV, releasing huge amount of entropy. By its own, the entropy release in 
the decay of flaton is not enough to resolve moduli problem. But the decay of LHu 
condensation before the end of thermal inflation provides an additional dilution. 
As the result, moduli can be dilute out to a safe level. Interestingly axion coupling 
constant and flaton decay temperature are tightly constrained to be fa ~ 10^^ GeV 
and Td ~ 100 MeV, respectively. Remarkably, model parameters in this case are 
likely to give naturally right amount of baryon asymmetry at present. 

The main component of dark matter in our model is expected to be the axino 
whose mass is nearly fixed to be about 1 GeV to match cold dark matter abundance 
at present. Axions from misalignment and strings also contribute to cold dark 
matter, but they are sub-dominant. Axions from the decay of flaton are sub- 
dominant too, but they are expected to be hot and highly relativistic at present in 
our model. The current number density spectrum of the hot axions may provide 
experimental and/or observational tests of our model in the future. The current 
number density spectrum of the hot axions may provide an observational test of 
our model in the future. 
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A Renormalization group equations 

The relevant part of superpotential of Eq. (1) for the running of soft SUSY- 
breaking parameters associated with (j) and is 

VFpart = XuLHuiy + ^A$$l/2 (10^) 

Then, one finds renormalization group equations (RGEs) 

= ^|A,PK + ™S + |A.P) (102) 
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dlnQ 


87r2 ' 


da^ 


1 


dlnQ 


167r2 


dX^ 


1 


dlnQ 


167r2 



(103) 

15 \ 

ya$|A$p + 4a^A$A* + 2a$|A^n (104) 

A$ (^|A$|2 + 2|A^|2^ (105) 

where Q is the renormahzation scale, mf is the soft mass-squared of a scalar field 
i, and Oj = AiXi with Ai the ^-parameter of the trilinear coupling associated with 
Aj. From Eq. (105), Eq. (104) can be rewritten in terms of \Aci,\ as 



d\Ais,\ 



\ [5|A4.|2|A<j|2 + 2 {A^Al + AlA,) |A,|2] (106) 



dlnQ 87r2 

If A$ ^> A,y, which should be the case for consistency of theory as described 
in the text, we can ignore all the contributions with \y. Then, Eqs. (102), (103), 
(106) and (105) can be approximated to 

^ = ^\X.mml + 2ml + \A,\^) (108) 
'^l^*'' - ^ 'A^lV-fl' (109) 



dlnQ 87r2' 

|A$P (110) 



dlnX^ 5,2 



dlnQ 327r2 

and in order for only to be negative at low renormahzation scale we need 

m2<m2 + |^^|2 , m2>|^$|2 (m) 

For example, for a potential only with the soft-mass term of (j), a numerical analysis 
shows that input values 

m^(Mpi) = 250 GeV , \A^{Mpi)\ = 150 GeV , m^{Mpi) = 100 GeV 

(112) 

with |A$(Mpi)| = 1 generates vacuum at (/iq ^ 6 x 10"'^2 GeV where 



m^~210GeV , |^$| ~ llOGeV , ^-ml = l2GeY (113) 
with \X^\ ~ 0.84 and mpq ~ 23 GeV the physical mass of (p. 
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